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ABSTRACT 



This work presents three dimensional, magneto-hydrodynamical simulations 
of the formation and early evolution of a subgroup of planetary nebulae 
that exhibit a variety of point-symmetric structures. For bipolar nebulae, 
the formation of point-symmetric nebular shapes along the inner borders of 
their opposing lobes and/or collimated outflows or jets internal or external 
to their cavities, is reproduced by considering a steady misalignment of the 
magnetic collimation axis with respect to the symmetry axis of the bipolar 
wind outflow, defined perpendicular to the equatorial density enhancement. 
Elliptical planetary nebulae with ansae displaced from the symmetry axis in 
point-symmetric fashion are reproduced through the same process by reducing 
the equatorial density enhancement. This mechanism represents an alternative 
explanation to some cases where morphological appearances give the impression 
of the action of a symmetric, rotating or precessing jet from the central source. 
The computational survey reveals that jet formation is detected only for dense 
enough winds with mass loss rates ;^ 10^^ M0 yr^"*^ . For lower mass loss rates 
the jets tend to vanish leaving behind only ansae-like structures at the tips of 
the lobes, as observed in some cases. The results are rather independent of the 
wind terminal velocity, since magnetized bubbles behave adiabatically for low 
wind velocities (~ 100 km s~^ ), which in the absence of a magnetic field would 
behave as momentum driven. 

Subject headings: Hydrodynamics — ISM: jets and outflows — planetary nebulae: 
general 
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1. Introduction 

Numerical simulations that reproduce the shapes of planetary nebulae (PNe) have 
helped to understand or give insights into the basic physics acting behind the formation 
of the different morphologies that are commonly observed in PNe. The two-wind model 
proposed by Kwok, Purton & Fitzgerlad (1978) has been applied and refined in a number of 
numerical simulations based on hydrodynamic models (e.g. Khan & West 1985; Mellema, 
Eulderink & Icke 1991; Frank & Mellema 1994; Dwarkadas, Chevalier & Blondin 1996). 
During the last decade, though, it has become evident that the production of highly 
collimated, bipolar, supersonic outflows, play a major role in the evolution of a substantial 
number, if not all, of PNe (e.g. Lopez 1997 for a review). 

The main nebular morphologies are usually well reproduced in the hydrodynamic 
models, however, the production of jets and ansae have met with certain difficulties; in 
particular, fast winds do not tend to converge into stable structures to form ansae and jets 
(e.g. Dwarkadas & Balick 1998). 

On a somewhat different approach. Chevalier & Luo (1994) have explored the effects of 
a rotating star with a magnetized wind on the formation of aspherical bubbles. Following 
this scheme, Rozyczka & Franco (1996) and Garcia-Segura et al. (1997a, 1999), performed 
2-D MHD simulations of this magnetized wind in cylindrical and spherical calculations, 
respectively, showing that magnetic tension may indeed be responsible for the generation of 
jets in PNe. Recently, Garcia-Segura (1997) (paper I) presented full 3-D MHD models of 
PNe where jets and ansae are convincingly reproduced as the result of magnetic coUimation 
of the post-AGB wind. 

A particularly intriguing case in PNe morphologies are those that display point- 
symmetric structures. Point-symmetry has been recognized in a wide variety of PNe 
(e.g. Guerrero, Vazquez & Lopez 1998). In some instances, where reliable kinematical 
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information is available, the presence of bipolar, rotating, episodic jets or BRETS have 
been inferred to explain similar structures (c.f. Lopez 1997). 

Cliffe et al. (1995) performed 3-D numerical simulations introducing in their 
computational mesh the presence of a BRET and reproducing the main morphological 
features observed, in particular for the case of Fleming 1 (Lopez, Meaburn & Palmer 1993). 
They also noted that some bipolar nebulae, such as Hb 5 present distinct, clumpy, S-shaped 
intensity distributions in their opposing lobes, which led them to argue that these structures 
could also be originated by a precessing jet. The numerical simulations presented in this 
paper have been mainly motivated by this group of objects, though the results expand 
a wide range of configurations, including the conditions for jet detection in PNe, ansae 
formation in point-symmetric modes and 'regular' bipolar and elliptical PNe presenting 
some degree of asymmetry. 

In § 2 we discuss the initial considerations and the motivation for the present paper. 
§ 3 describes the methodology. § 4 presents results for bipolar nebulae and § 5 for ellipticals. 
§ 6 Discusses the results and a summary is given in § 7. 

2. Initial considerations 

As an illustration of the main type of morphological structures that are the focus 
of the present simulations, a mosaic of images of these type of nebulae is presented in 
Figure 1. The objects included in this representative sample are, (la) Hb 5, (lb) He 2-103, 
(Ic) He 2-429 and (Id) NGC 2371. Hubble 5 is a bipolar nebula with distinct internal 
S-shaped intensity distribution. Recent HST images of this object obtained by Balick, Icke 
and Mellema show a great deal of detail in the secondary, tilted outflow. He 2-103 is an 
elliptical nebulae where the equatorial density enhancement is clearly not perpendicular to 
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the elliptical shell; furthermore, faint nebular extensions are discerned displaced from the 
main elliptical axis on both sides of the nebula. He 2-429 (see Guerrero et al. 1998 and the 
lAC Catalog, Manchado et al. 1996) is an elliptical nebula with twisted, point-symmetric 
extensions protruding outside the main nebular shell. Finally, NGC 2371 has a dumb-bell 
appearance with intensity maxima in the equator defining a minor axis which is not 
orthogonal to the main nebular axis. 

The main, standard, working hypothesis that sets the initial conditions of the 
simulations is that the equatorial density enhancement (EDE) produced in the late AGB 
phase defines the main symmetry axis of the bipolar/elliptical wind outflow, i.e. the EDE 
is orthogonal to the main bipolar /elliptical structure. 

In paper I, in order to reproduce point-symmetric structures, a wide binary nucleus 
was considered and the magnetic coUimation axis (i.e. the rotation axis of the primary 
star) was assumed to precess, with a period of the same order of the kinematic age of the 
nebulae (~ lO'^ yr). 

In the present paper, instead of invoking an active rotation of the magnetically 
collimated outflow, a misaligned conflguration is produced by introducing a steady tilt of 
the magnetic collimation axis with respect to the bipolar /elliptical wind outflow. Hence, 
this is a particular case of paper I, in which the precession period is nule or much larger 
than the kinematic age of the nebula. 

Thus, in order to achieve this steady misalignment, the rotating core of the AGB, 
which will provide the magnetized fast wind, must have a rotation axis tilted with 
respect to the axis deflned by the EDE. Our current limited understanding of the detailed 
processes that produce the EDE at the AGB phase precludes a full physical justiflcation 
for this assumption. Binary nuclei could represent the best way to attain this effect given 
the diversity of paths that they can take during their evolution through the AGB (e.g. 
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Yungelson, Tutukov & Livio 1993). For example, wide binary systems may induce this 
misalignment since the ejected envelope's rotation axis of the primary star could be tidally 
altered by the secondary and the EDE formed in the orbital plane, with a resulting tilt with 
respect to the rotation axis of the primary. Alternatively, systems which undergo a common 
envelope evolution, in which the rotation axis of the primary star is not orthogonal to the 
orbital plane could also produce the same effect. Close binary systems with substellar 
companions that may lead to the formation of accretion disks (Reyes-Ruiz & Lopez 1999) 
and mergers could also play a significant role here. 

However, beyond these arguments, which badly need additional investigation, the best 
initial justification that has motivated at this time further exploring of this assumption is 
the apparent asymmetry observed in the equatorial regions of a substantial number of PNe 
and the corresponding point-symmetric condensations and shapes (see Figure 1). 

The numerical models do not consider an explicit physical form for the source of the 
wind, just an outflowing magnetized wind from a spherical radius of r = 0.005 pc after t=0. 
Regions below this limit lie beyond the current resolution of the simulations. 

The general idea behind these MHD models is sketched in Figure 2. Here, a rotating 
star produces the outflowing magnetized wind, in which the toroidal component of the 
magnetic field is dominant respect to the poloidal contribution at large distances. 

3. Three-dimensional MHD models 

The simulations have been performed using the magneto-hydrodynamic code ZEUS-3D 
(version 3.4), developed by M. L. Norman and the Laboratory for Computational 
Astrophysics. This is a finite-difference, fully explicit, Eulerian code descended from the 
code described in Stone & Norman (1992). A method of characteristics is used to compute 
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magnetic fields as described in Clarke (1996). ZEUS-3D does not include radiation transfer, 
but we have implemented a simple approximation to derive the location of the ionization 
front for arbitrary density distributions (see Bodenheimer et al. 1979; Franco et al. 1989, 
1990, Garcia-Segura & Franco 1996). This is done by assuming that ionization equilibrium 
holds at all times, and that the gas is fully ionized inside the H ll region. 

The location of the ionization front in any given direction (^, 0) from the photoionizing 
source is given by J n^{r,6,(f))r'^dr ^ F^/AiraB, where symbols have their usual meaning. 
The models include the Raymond & Smith (1977) cooling curve above 10"^ K. For 
temperatures below 10*^ K, the shocked gas region is allowed to cool down with the radiative 
cooling curves given by Dalgarno & McCray (1972) and MacDonald & Bailey (1981). 
Finally, the photoionized gas is always kept at 10"^ K, so no cooling curve is applied to 
the HII regions (unless, of course, there is a shock inside the photoionized region). The 
minimum temperature allowed in all models is set to 10^ K. The computations are done in 
Cartesian coordinates (x, y, z), with outflowing outer boundary conditions. 

The models have grids of 100 x 100 x 200 equidistant zones, with an extent of 
0.1 X 0.1 X 0.2 pc in x, y and z respectively. The wind, whose origin is centered at the 
center of the computational volume, is set based on the rotating wind solutions given 
by Bjorkman & Cassinelli (1993), Ignace et al. (1996, 1998) and the approach to those 
equation given by Garcia-Segura et al.(1999). These equations permit the introduction of 
the wind self-consistently. For the reasons given in the preceding section, in the present 
calculations we do not address a particular functional shapes for the EDE, but simply state 
its presence independently of its origin. 

The magnetic field of an outflowing wind from a rotating star can be described by 
two components, B^ and Bj- (Chevalier & Luo 1994; Rozyczka & Franco 1996; paper I; 
Garcia-Segura et al. 1999). The radial fleld component can be neglected since Bj. ~ r^^ 



while the toroidal component -B^ ~ r~^, and the field configuration obeys V ■ B = . In all 
the simulations, the toroidal field is set always perpendicular to the z axis, therefore only 
B^^ and B^ are introduced in the source wind sphere, with B^^ = 0. 

As in paper I, we will use two useful dimensionless parameters containing the basic 
information of the models. The first one Q is the velocity ratio of the stellar rotation to the 
critical rotation, and the second parameter a is the ratio of the magnetic field energy density 
to the kinetic energy density in the fast wind. The wind compression field ("WCField") 
found by Ignace et al.(1998) is also introduced. However, since Q is less than 10% for the 
fast wind in all the calculations, this correction is of little relevance in these cases. 

The AGB wind is set with a Voo = 10 km s~^ and a M = IO^^Mq yr^^ as a standard 
value in all the models. The AGB wind is set unmagnetized in all models, because the 
estimated values for a are very low in this regime (see paper I). Since we are interested in 
the early phases of the PN evolution, the velocity of the fast wind is set to 100 km s^^ . 
This value should reasonably match the conditions at the departure from the AGB phase. 

A more realistic approach to this problem would need to include the evolving wind 
conditions with gradual transitions of the wind from v^o ~ 10 km s~^ up to ~ 1000 
km s~^ (or higher), and mass-loss rates varying from M = 10^® Mq yr~^ (or higher) down 
to 10"^ (or lower). In this case time dependent effects invoke additional behaviors on the 
dynamics (e.g. Dwarkadas & Balick 1998). They are out of the scope of the present paper 
and are deferred to a future work. 

The parameters of the runs are summarized in Table 1. These include, an unmagnetized 
reference case (model AO); bipolars with different inclinations of the bipolar wind outfiow 
symmetry axis with respect to the magnetic coUimation axis and two different mass loss 
rates (models B and C, respectively) and ellipticals with a single constant tilt and three 
different mas loss rates (models E) 
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4. Bipolar Nebulae 

Bipolar nebulae are mainly characterized by an equatorial waist with a substantial 
pole to equator density contrast. In order to achieve such a density contrast, we have used 
a value of f] = 0.98 in the set up, which is obtained with a rotational velocity of 6.929 
km s~^ . 



4.1. The unmagnetized bipolar case 

As a starting point, the response of a spherically symmetric, unmagnetized wind 
starting from a spherical shape is shown in model AO (Fig. 3). This unmagnetized case 
shows how the wind behaves in such a medium, i.e., the wind is impeded to expand in the 
equator due to the EDE, but it easily expands along the polar direction. This would be the 
classical hydrodynamic case that has been thoroughly studied in the past (c.f. see Frank 
1999 for a review on bipolar outflows). Because of the radiative conditions of the reverse or 
terminal shock of such a wind in this set up ( tpost-shock ~ 10^ K for f shock = 100 km s~^ ), 
the bipolar bubble is momentum driven. This particular dynamics generates ram-ram 
pressure instabilities ( Garcia- Segura et al. 1996; Garcia-Segura et al. 1997b), however, 
given the low resolution of the grid, these instabilities are not reproduced in the present 
study, and only some indications of such a behavior is present as the small corrugation that 
appears in the lobes of model AO. 

Model AO is important as a reference model since it shows that fast unmagnetized 
winds, starting from spherical conditions and expanding within an EDE cannot converge to 
form jets or ansae. In more adiabatic conditions, the lobes become two tangential spheres 
(Garcia-Segura et al. 1999). It is precisely this behavior that rules under the presence of an 
EDE, i.e. the bipolar lobes will grow always perpendicular to the EDE, and is the EDE the 
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responsible for the equatorial waist observed in bipolar nebulae. 

4.2. Magnetized Bipolar Nebulae 

Keeping in mind the behavior of the reference model AO, a magnetized wind is now 
introduced to visualize the strong change that effects the dynamics. Model BO (Fig. 3) 
only differs from model AO in the magnetized wind, with a o" = 0.01. But the result is 
remarkable, a pair of jets form at both sides in the lobes. The explanation of this behavior 
is fully described in Rozyczka & Franco (1996), and it results mainly from the action of 
the hoop stress caused by the toroidal magnetic field. The rest of the bubble behaves now 
adiabatically because the magnetic pressure does not allow the gas to pile up towards 
the inner edges of the lobes. This is also the main reason why model BO expands much 
faster than model AO, since energy conserving bubbles in power laws of the type r~^ have 
different dynamics than momentum conserving bubbles (Garcia-Segura & Mac Low 1995 
and references therein). 

Although the hoop stress is always present in this type of flows, not always the 
conditions are in favor to form observable jets. Such a case is computed in model CO 
(Fig. 3), in which the mass- loss rate of the fast wind has been decreased by an order of 
magnitude. In such conditions, the fast wind piles up at the poles, but the radiative cooling, 
proportional to n^, is not sufficiently strong to allow the formation of visible jets, and 
ansae-like structures are formed instead. This difference is more evident in the synthetic 
pictures shown in figure 7 for models Bs and Cs, which are discussed further below. 

As mentioned before, the symmetry axis of the bipolar wind outfiow is considered 
orthogonal to the EDE. Models BO and CO are computed with the magnetic collimation 
axis also orthogonal to the EDE. The following models break this condition by introducing 
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different tilts between these axes to explore the resultant effects, which are shown to lead 
to point-symmetric shapes. 

The tilt is introduced at the EDE for simplicity. Thus, the magnetic collimation 
axis remains parallel with the z axis, and the same set up described above is used for the 
magnetized wind. 

In order to have a general picture of the effect caused by the tilt, we have computed 
several models in which the tilt angle varies from low values, 5° (models B5 and C5), 
middle values, 15° (models B15 and C15) and large values, 45° (models B45 and C45). The 
evolution of the two last models is shown in figures 4 and 5, respectively. 

The dynamics of the bubbles in such conditions are quite complex and the limited 
resolution of the models hinder a full, quantitative description. Nevertheless, the effect of 
the tilt can be summarize qualitatively as follows. 

At the beginning of every simulation, the spherically symmetric shocked flow is driven 
by the magnetic pressure gradient after crossing the reverse or terminal shock and is 
launched in the z direction on both sides of the wind source region. Since the magnetic 
field is toroidal, the magnetic pressure is nule towards the z axis in the free expanding wind 
(but not in the shocked region), while it is maximum in the equatorial plane, consequently 
building a large pressure gradient. Gradually, the flow notices the effect of the EDE and 
is deflected perpendicular to it. The magnetically pressurized bubble is distorted, and the 
flow, which now follows the tendency of the bipolar cavity, pushes and deflects the column 
of gas in which the hoop stress is larger. When this column is dense enough, the formed 
jets are deflected an acquire the appearance of precessing jets, which eventually will collide 
with the bipolar cavity (see e,g, model B45). Thus, there are two competing effects, one is 
the tendency of the gas to flow perpendicular to the EDE, and the other is the tendency of 
the gas to flow parallel to the z axis as imposed by the magnetic pressure and the hoop 
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stress. When the tilt is nule, i.e., when both effects do not compete with each other, the 
bubble acquires the largest expansion velocity at the polar regions. 

5. Magnetized Elliptical Nebulae 

Elliptical nebulae are characterized by a moderate pole to equator density contrast. 
We have adopted a value oi Q = 0.80 for the AGB wind, which is obtained with a 
rotational velocity of 5.65 km s^^ . The AGB wind is set with a foo = 10 km s~^ and a 
M = 10~^Mq yr"-*^ as a standard value for models El and E2, while M = 1O~^M0 yr~^ 
is used for E3 and E4. The mass loss rate of the fast wind is varied from values of 
M = 10-^Mq yr~^ (El) to M = IO^^Mq yr'^ (E3 and E4). A tilt of 15° is used in this 
case for all models for simplicity. 

An example of the model evolution is shown in Fig. 6 for model El. The evolution 
is similar to the models described above for bipolar nebulae, with the difference that 
the bubble grows further at the equator since the EDE is not as restrictive as before. 
Figure 7 (right panels) shows the apparent differences between models El, E2, E3 and 
E4, dictated by varying the mass loss rates of the fast wind from relatively large to lower 
values. The polar, piled-up gas in model El (M = 1O~'^M0 yr^^ ) is able to cool down 
efficiently in a dynamical time scale comparable to the computed time and the axial jets 
become readily apparent. In model E2 (M = IO^^Mq yr^^ ) the cooling is less efficient 
but still present, giving rise to the formation of ansae-type structures. Models E3 and E4 
(M = IO^^Mq yr^^ ) do not show any apparent signs of piled-up gas (see Figure 7). It 
should be noted, however, that given the limited numerical resolution of these models, 
the cooling regions are not resolved and the mass loss rates should be taken as estimative 
values. However, since these models have energy conserving dynamics (as imposed by the 
magnetic pressure) rather than momentum conserving, these estimations for the mass loss 
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rates are quite independent of the terminal wind velocities. In the absence of a magnetic 
field, winds of the order of 100 km s~^ have strongly radiative terminal or reverse shocks 
(see for example Frank et al. 1996), and the radiative conditions are strongly dependent on 
the wind velocity and not only of the mass loss rate. 

6. Analysis of the models 

In order to test the numerical results of the models, it would be convenient to 
run convergence tests at lower and higher resolution. The latter becomes prohibitively 
demanding for our computational resources. An experiment at half the resolution, which 
gives a qualitative indication on the reliability of the model is presented in Figure 6b. This 
figure shows a comparison for model El with two different resolutions. Althought, the 
individual details in both simulations in that figure differ by obvious reasons, the main 
features, such as the shape of the bubble and the formation and extension of the jets remain 
rather similar. The densities are larger in the high resolution run, indicating that in order to 
make quantitative estimates of physical parameters, one needs twice the resolution at least 
to be sure that the dynamics has been resolved, i.e., to find a convergence. For example, 
what is the post-shock cooling length and how does this compare with the resolution of the 
simulations ? 

In order to resolve a cooling distance, one needs basically to get rfcooi > x^-aax/n^ , 
where x^ax is the physical size of the grid, and n^ the number of zones (i.e. rfcooi bigger 
than one zone). However, several zones are needed to resolve detailed cooling processes 
(e.g., any sort of microphysics). A rough estimate of the cooling timescales, can be obtained 
from equilibrium line cooling. In this case, the cooling time of a parcel of shocked gas is 
given by tcooi = 3^/^ Vg/{2^ A po) — 0.3 t;| Hq^ yr (Franco et al. 1994), where vs = "Ws/lO^ 
cm s~^ is the shock velocity, po is the pre-shock mass gas density and tiq is this density 
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in units of 10® cm^^, A = 3.9 x 10^^ Qp^^'^ cm® s^^ g^^, and Qp is the mass per particle in 
units of rriH {i- e., Qp = 1.4/2.3 for a fully ionized gas with cosmic abundances). For typical 
numbers in our simulations {vs ~ 100 km s~^ , n ~ 10 cm^^) we obtain tcooi ~ 30 yr, which 
give us (icooi ~ 2 X 10^^ cm, just slightly smaller than the zone zise of 3 x 10^^ cm. One 
should also consider that in MHD models conditions are such that cooling distances are 
more favourable resolved perpendicular to the field lines given the tendency to adiabatic 
behaviors in this case. 

Since the MHD simulations behave similar to the hydro solutions along the Z axis 
(parallel to the field lines), except for the extra pressure (hoop stress) that the magnetic 
field produces towards the rotation axis of the central star, it is possible to use the relations 
of conical converging flows (Canto, Tenorio-Tagle & Rozyczka 1988; or the specific solutions 
for PPNe by Frank, Balick & Livio 1996), as an approximation for the parameters in the 
jets. This is in fact a conical converging, magnetized flow problem. Thus, for example, 
the density in the jets can be estimated by following formula (23) in Canto, Tenorio-Tagle 
& Rozyczka (1988), where pj ~ 6 x lO^sm"^^ [fo/100 km s~^ ] po • For the case of our El 
model, one gets (for 6 ~ 50°) pj ~ 1 x lO^po- On the other hand, the computed models give 
Pj ~ 1 X lO^po for the low resolution run and pj ~ 1 x lO^po for the higher resolution one. 
Since the resolution is not big enough (as discused above), and in both cases, po ~ 1 x 10^^^ 
g cm^^ , it is expected to find a jump of three orders of magnitude on the density ratio at 
twice the resolution, which would match the analytical value. 

7. Discussion 

The computational survey reveals that the formation of jets and ansae, or none, can be 
achieved by varying the mass-loss rate within a common framework of magnetized winds. 
In principle this would allow to have an estimate of the mass-loss rates involved in the 
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formation of the different types of collimated outflows found in PNe. In particular, it is 
found that for the formation and detection of jets, the mass- loss rates must have values 
M ^ 10"'^ M0 yr^-*^ . In the case of ansae a M ~ 10~^ M0 yr~^ is required. Below these 
values, i.e. for typical cases of M ~ 10~^ M© yr~^ or lower, no collimated features are 
formed or revealed by the models. These numbers, however, must be taken with caution 
given the low resolution used in these 3-D calculations. Nevertheless, there are additional 
studies that indicate that jets in PNe must be formed during a brief lapse of high mass-loss 
rate in the pre-planetary nebula stage. For example, Reyes- Ruiz & Lopez (1999) have 
concluded that when an accretion disk is involved in the formation of collimated outflows 
in PNe, the necessary high accretion rates, M ~ 10"^, needed to launch an axis-symmetric 
wind (as in the case of young stellar objects, e.g. Calvet 1998) are attained only for a few 
thousand years during the late AGB or pre-planetary phase. 

Figure 7 shows different models in which jets have formed. The projected shapes in 
the tilted cases give the impression of a varying direction in the outflow, closely resembling 
the point-symmetric structures often observed in a number of PNe. However here, no 
active precession or rotation in the source of the outflows is involved. These shapes result 
by the steady misalignment of the magnetized wind's axis with respect to the axis of the 
bipolar wind outflow. This configuration can lead to a hydrodynamical deflection of the 
magnetized wind on the bipolar cavity. This is a novel feature revealed by the simulations 
and which can be tested in the future by detailed kinematical and spectropolarimetric 
studies. For example, in this case the deflected gas is expected to show a kinematic behavior 
that deviates from the rest of the nebular expansion and line emission in the magnetically 
collimated regions to show polarization effects. Hydrodynamical deflection can only account 
for S-shaped morphologies, helicoidal type structures would require some sort of rotation or 
wobbling from the wind source. 
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Although the PNe that were originally classified as point-symmetric comprise a very 
reduced group, the presence of point-symmetric structures in PNe has been shown now to 
span a much wider sample (Guerrero, Vazquez & Lopez 1998). For example, signatures 
of point-symmetry can be appreciated in nearly 20% of the sample contained in the lAC 
morphological catalog (Manchado et al. 1996). Furthermore, practically every PN that has 
been recently observed at high spatial resolution with the Hubble Space Telescope shows 
some degree of point-symmetry, indicating that the mechanisms that produce it should be 
rather common. The representative models shown in Figures 2-7 reveal the variety of 
point-symmetric configurations that can be formed under the assumptions and conditions 
of the present study. The inclusion of binary nuclei in this analysis can only but strengthen 
the conditions for the development of point-symmetry. 

8. Summary 

We have made a computational survey of 3-D MHD simulations of young planetary 
nebulae in which are introduced different steady misalignments of the magnetic collimation 
axis with respect to the symmetry axis of the bipolar/elliptical wind outfiow. The latter 
defined as perpendicular to the equatorial density enhancement. This steady tilt can 
also be interpreted as the result of a very large precession period. The simulations show 
that in these cases a hydrodynamical defiection of the magnetized collimated wind on the 
bipolar/elliptical cavity can produce morphologies which may resemble the presence of 
precessing or rotating sources. In this way, the resulting morphologies are able to reproduce 
a considerable number of PNe which present some degree of point-symmetry in their 
structures. The survey also suggests that the formation and detection of jets and ansae, or 
none, strongly depends on the mass-loss rate of the driving, magnetized wind. In particular, 
the presence of jets is detected in the models only for mass-loss rates M ^ 10"^ M© yr~^ ; 
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below these rates the jets tend to vanish in the models, keeping for M ~ 10^^ M© yr^^ 
only ansae type structures. 
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Figure Captions 

Figure 1. A representative sample of the type of nebulae that are the focus of this study. 
1(a) Hb 5; 1(b) He 2-103; 1(c) He 2-429; 1(d) NGC 2371. See models B45, E3, El and 
E4 in Figure 7, respectively for comparison. 

Figure 2. Schematic view of the models. The magnetic collimation axis (dashed line) 
is tilted with respect to the symmetry axis of the bipolar wind outflow, which is 
perpendicular to the equatorial density enhancement. 

Figure 3. Evolution of the logarithm density distribution (g cm~'^) of a two-dimensional 
slice in the x-z plane which crosses through the middle of the computational volume. 
Top panels correspond to the unmagnetized case, model AO (M = 10^^ M0 yr^^ ) 
for t = 1250 and 2500 yr. Middle panels correspond to the magnetized case, model 
BO {a = 0.01 and M = IO^^Mq yr^^ ) for t = 900 and 1800 yr. Bottom panels, 
magnetized case, model CO {a = 0.01 and M = 10~^ Mq yr~^ ) for t = 1250 and 
2500 yr. In all these cases the symmetry axis of the bipolar wind outflow and the 
magnetic collimation axis are coincident. 

Figure 4. Evolution of the logarithm density distribution (g cm^^) of a two-dimensional 
slice in the x-z plane which crosses through the middle of the computational volume 
for model B45 (o" = 0.01 and M = 10^^ Mq yr^^ ). Here the symmetry axis of the 
bipolar wind outflow is tilted by 45° with respect to the magnetic collimation axis (t 
= 350, 700, 1000, 1350, 1650, 2000 yr). 

Figure 5. Same as Figure 4 but for model C45 where M = 10~^ Mq jr~^ (t = 400, 850, 
1250, 1650, 2100, 2500 yr). 

Figure 6. Same as Figure 4 for model El {a = 0.01 and M = 10~^ Mq yr^^ ). Here 
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the symmetry axis of the bipolar wind outflow is tilted by 15° with respect to the 
magnetic coUimation axis (t = 250, 550, 800, 1050, 1350, 1600 yr). 

Figure 6b. Two-dimensional slices of Model El for 50 x 50 x 100 (top) and for 
100 X 100 X 200 (bottom) equidistant zones in x, y and z respectively. 

Figure 7. Two-color, composite, synthetic images of the models where red represents the 
shock-dominated (non-photoionized) material and green the photoionized gas. The 
line of sight corresponds to the y axis, being the x-z plane the plane of the sky. Bipolar 
models in the left panels show the magnetized case with no tilt between the symmetry 
axis of the bipolar wind outflow and the magnetic collimation axis (model BO) and 
the magnetized cases for three different tilts (B5, B15 and B45). In these models 
M = 10^^ Mq yr^^ and a = 0.01. The central panels correspond to the bipolar 
models, as before, but with M = 10^^ M0 yr^^ . Elliptical models are shown in the 
right panels. A tilt of 15° is present in these four models. El (M = 10"'^ M0 yr"-*^ , 
a = 0.01), E2 (M = 10"^ M© yr'^ , a = 0.01), E3 (M = 10^^ M© yr'^ , a = 0.01) 
and E4 (M = 10^^ Mq yr~^ , a = 0.99). Each run in the figure corresponds to: BO 
and B5, 1800 yr; B15 and B45, 2000 yr; CO, C5, C15 and C45, 2500 yr; El, 1600 yr; 
E2, 1550 yr; E3, 2000 yr and E4, 350 yr. 
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Table 1: Wind Properties 
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